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5t 1 ference patterns are produced without making reference to the solution of a wave equation 
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1. Introduction 

In 1802, Young performed a double-slit experiment with light in order to resolve the ques- 
tion whether light was composed of particles, confirming Newton's particle picture of light, 
or rather consisted of waves. ^' His experiment showed that the light emerging from the slits 
produces a fringe pattern on the screen that is characteristic for interference, discrediting 
Newton's corpuscular theory of light. -'^-* Hence, from the point of view of classical physics, the 
particle and wave character of light did not seem to be compatible. Moreover, the interpreta- 
tion in terms of particles or waves of the observations in experiments with light became even 
more complicated after conduction of the Michelson-Morley experiment^) which provided evi- 
dence that light waves do not need a medium (the ether) to propagate through, in contrast to 
water and sound waves which require media. However, explanation of the photoelectric effect 
by Einstein in terms of photons^' is perhaps the most direct and convincing evidence of the 
corpuscular nature of light. Einstein's explanation of the photoelectric effect was the start of 
understanding the quantum nature of light and influenced the development of the concept of 
wave-particle duality in quantum theory. 

In 1924, de Broglie introduced the idea that also matter, not just light, can exhibit wave- 
like properties.^-* This idea has been conflrmed in various double-slit experiments with mas- 
sive objects such as electrons, ^"^^ neutrons,^' ^'^^ atoms^^'^^^ and molecules such as Cqq and 
Cyo/^'^^) all showing interference. In some of the double-slit experiments^' '^'^^•* the inter- 
ference pattern is built up by recording individual clicks of the detectors. Identifying the 
registration of a detector click, the "event" , with the arrival of a particle and assuming that 
the time between successive clicks is sufficiently long such that these particles do not interact, 
it becomes a challenge to explain how the detection of individual objects that do not interact 
with each other can give rise to the interference patterns that are being observed. Accord- 
ing to Feynman, the observation that the interference patterns are built up event-by-event is 
"impossible, absolutely impossible to explain in any classical way and has in it the heart of 
quantum mechanics". ^^' 

Although wave-particle duality is a central concept of quantum theory, in practice quantum 
theory only works with wave functions to describe the total system under study. In order to 
describe the single occurrences observed in various experiments the process of wave function 
collapse has been introduced. However, the precise mechanism of a wave function collapse is 
not known. 

Recently, various experiments have been performed that measure individual events gener- 
ated by microscopic objects. Hence, it is of interest to study how the particle and wave picture 
of these experiments are contradicting each other. It is often said that wave properties like 
interference cannot be realized by non-interacting particles which satisfy Einstein's criterion 
of local causality. In earlier work we have presented an event-based corpuscular simulation 
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model which demonstrates that such particles can indeed produce interference patterns and 
applied it to a variety of single-photon experiments like beam splitter and Mach-Zehnder in- 
terferometer experiments, Wheeler's delayed choice experiments and many others. ^^"^^•* What 
these experiments have in common is that the interference can be described as two-path inter- 
ference, that is the observed interference pattern is the result of having only two possible paths 
for the particles travelling between source and detector. In order to simulate such experiments 
it is sufficient to use adaptive models for the optical apparatuses and to use detectors that 
simply count the number of detection events. ^^"^^-* In this paper we extend the simulation 
model towards simulating multipath interference patterns as observed in single-photon two- 
beam interference and two-slit experiments, for example. Detectors that are simply counting 
the detection events cannot be used for this purpose. Therefore we introduce a new simulation 
model for the single-photon detector that takes into account the memory and threshold behav- 
ior of such a detector. The model is a natural extension of the earlier work mentioned and is 
fully compatible, that is interchanging in our earlier work the simple counting detector model 
with this more complex detector does not change the conclusions. In this sense, the present 
detector model adds a new, fully compatible, component to the collection of event-by-event 
simulation algorithms. 

Note that the event-based simulation model is not a corpuscular model in the classical- 
physics sense. In our model, particles are objects that carry information. As a particle en- 
counters a material device, it exchanges information with this device. In our model, this in- 
formation exchange is the cause of the appearance of an interference pattern. In other words, 
in our approach we construct a mechanism which produces wave-like phenomena by local 
variables only. To this end, we introduce independent objects which carry information. These 
objects we call "particles". Each particle interacts with the material of the device only and 
the effect of many of such interactions is to build up a situation which causes the appearance 
of a first-order interference pattern. 

To head off possible misunderstandings, the present paper is not concerned with an inter- 
pretation or an extension of quantum theory nor does it affect the validity and applicability 
of quantum theory. Furthermore, the event-based detector models that we introduce in this 
paper should not be regarded as realistic models for say, a photomultiplier or a photographic 
plate and the chemical process that renders the image. Our aim is to show that, in the spirit 
of Occam's razor, these very simple event-based models can produce interference patterns 
without making reference to the solution of a wave equation. 

Although waves can be the physical cause of interference, the key point of our work is 
that it is wrong to think that waves are the only possible physical cause of interference: In 
our approach, the clicks produced by non-interacting / non-communicating detectors, caused 
by non-interacting / non-communicating particles that arrive at single detectors one-by-one. 
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Fig. 1. Schematic diagram of an interference experiment with a Fresnel biprism (FBP).^'^' S, Si, S2 
denote the point source and its two virtual images, respectively. The grey area is the region in 
which an interference pattern can be observed. 



build up a pattern that is identical to the one that is obtained by solving a wave equation. 
However, our event-based simulation approach does not require knowledge of the wave am- 
plitudes obtained by first solving the wave mechanical problem or requires the solution of the 
Schrodinger equation. Interference patterns appear as a result of an event-by-event simulation 
of classical, locally causal, adaptive dynamical systems. 

The paper is structured as follows. In Section 2, we introduce the interference experiments 
that we simulate. In Section 3, we review the main features of the photon detection process. 
Section 4 specifies the new detector models and the simulation algorithm in full detail. A 
Mathematica implementation of this algorithm for the case of the double-slit experiment can 
be downloaded from the Wolfram Demonstration Project web site.^^-* In Section 5, we compare 
the event-by-event simulation results with the numerical results obtained from wave theory 
for the two-beam interference experiments discussed in Section 2, showing that our event- 
based, particle-like approach reproduces the results of quantum theory without making use 
of concepts thereof. In Section 6, we propose a realizable experiment to test our event-based 
models for interference. Our conclusions are given in Section 7. 

2. Two-beam interference 

In this paper, we focus on interference experiments with single-photons, leaving the case 
of massive particles for further research. As a prototype problem, we consider two-beam 
interference experiments with a Fresnel biprism. "^'^' A schematic diagram of such an experiment 
is shown in Figure 1. A pencil of light, emitted by the source S, is divided by refraction into two 
pencils. ^''•' Interference can be obtained in the region where both pencils overlap, denoted by 
the grey area in Fig. 1. As a Fresnel biprism consists of two equal prisms with small refraction 
angle and as the angular aperture of the pencils is small, we may neglect aberrations.^*^' The 
system consisting of the source S and the Fresnel biprism can then be replaced by a system 
with two virtual sources Si and S2, see Fig. 1. Alternatively, following Young"^*^' we can let 
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Fig. 2. Schematic diagram of a sim.plified double-slit experim.ent with two sources 5*1 and ^2 of 
width a, separated by a center-to-center distance d, emitting light according to a uniform current 
distribution (see Eq. (1)) and with a uniform angular distribution, /? denoting the angle. The light 
is recorded by detectors D positioned on a semi-circle with radius X. The angular position of a 
detector is denoted by 6. 
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Fig. 3. Schematic diagram of a two-beam interference experiment with two line sources 5*1 and S2 
having a spatial Gaussian profile (sec Eq. (3)), emitting light according to a uniform angular 
distribution, j3 denoting the angle. The sources are separated by a centcr-to-centcr distance d. 
The light is detected by detectors D positioned at {X,y). 



the light impinge on a screen with two apertures and regard these apertures as the two virtual 
sources ^i and S2, see Figs. 2 and 3. Results of a single-photon interference experiment with 
a Fresnel biprism and a time-resolved interference experiment for the system schematically 
depicted in Fig. 3 are reported in Refs.^^' and,^^-* respectively. 

For all these simplified systems, a straightforward application of Maxwell's theory yields 
the intensity at the detection screen. We consider a few representative cases for which closed- 
form expressions can be obtained: 

• The sources Si and ^2 are lines of length a, separated by a center-to-center distance d, 
see Fig. 2. These sources emit light according to a uniform current distribution, that is 



J{x, y) = 5{x) [e(a/2 - \y - d/2\) + Q{a/2 -\y + d/2\)] , 



(1) 
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where G(.) denotes the unit step function. In the Fraunhofer regime, the hght intensity 
at the detector on a circular screen is given by^*^-* 

(qa sin 9 \ i ■ n 

sm ^ — \ r, qd sm 6 

_^j eos^^^— , (2) 

where ^ is a constant, q is the wave number, and 9 denotes the angular position of the 
detector D on the circular screen, see Fig. 2. 

• The sources 5*1 and 52 form a line source with a current distribution given by 

J{x,y) = 5ix)Y,e-^y-^'/'^'/'^\ (3) 

s=±l 

where a is the variance and d denotes the distance between the centers of the two sources, 
see Fig. 3. The intensity of the overlapping pencils at the detector reads 

I{y) = B (cosh ^ + cos (1Z|^^ ,-Ky^ +^V ^) / <^^ , (4) 

where i? is a constant, b = q^(T'^/{X'^ + q^cr'^)-, and (X, y) are the coordinates of the 
detector D (see Fig. 3). Closed- form expression Eq. (4) was obtained by assuming that 
d<^ X and a <^ X. 

• The two sources Si and 5*2 are circles with a radius a and their centers are separated by 
a distance d. The current distribution is given by 

J{x, y, z) = 6{x) [e(aV4 - {y - d/2f - z^) + e(aV4 - {y + d/2f - z^)] . (5) 

In the Fraunhofer regime, the light intensity at a detector placed on a sphere is given 

by30) 

iie) = c(^-^^^^^^^)\os''-^^^, (6) 

^ ^ V qasmO J 2 ' ^ ^ 

where C is a constant, 6 denotes the zenith of the detector D on the spherical detection 
screen and Ji(.) is the Bessel function of the first kind of order one. 

From Eqs. (2), (4) and (6), it directly follows that the intensity distribution on the detection 
screen, displays fringes that are characteristic for interference. 

3. Event-by-event simulation and detector model 

Imagine that individual particles build up the interference pattern one by one and exclude 
the possibility that there is direct communication between the particles (even if one particle 
has arrived at the detector while another particle is at the source or at a detector) . If we then 
simply look at Fig. 2 or 3, we arrive at the logically unescapable conclusion that the interference 
pattern can only be due to the internal operation of the detector: There is nothing else that 
can cause the interference pattern to appear. 

Obviously a simple, passive detector model that only counts the number of particles fails 
to reproduce the interference patterns of two-beam interference experiments in which there are 
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sources and detectors only, as in Figs. 2 and 3. Before we introduce new event-based models 
for the detector, it is expedient to review the conventional theory of the photon detection 
process. 

In its simplest form, a light detector consists of a material that can be ionized by light. 
The electric charges that result from the ionization process are then amplified, chemically in 
the case of a photographic plate or electronically in the case of photo diodes or photomulti- 
pliers. In the wave-mechanical picture, the interaction between the incident electric field E 
and the material takes the form P • E, where P is the polarization vector of the material. ^'^•' 
Treating this interaction in first-order perturbation theory, the detection probability reads 
Pdetectionit) = /g /(J((E^(t') • K(t' - t") ■ E{t")))dt' clt" where K(t' - t") is a memory kernel 
that is characteristic for the material only and ((.)) denotes the average with respect to the 
initial state of the electric field. ^^•' Both the constitutive equation^") P(^) = x('^)E(t<j) as 
well as the expression for Pdetectionit) show that the detection process involves some kind of 
memory. Furthermore, very sensitive photon detectors such as photomultipliers and avalanche 
diodes are trigger devices, meaning that the recorded signal depends on an intrinsic thresh- 
old. Conceptually, the chemical process that renders the image encoded in the photographic 
material plays a similar role. 

From these general considerations, it is clear that a minimal model for the detector should 
be able to account for the memory and the threshold behavior of the detectors. An event- 
based model for the detector cannot be "derived" from quantum theory, simply because quan- 
tum theory has nothing to say about individual events. '^^' Therefore, from the perspective of 
quantum theory, any model for the detector that operates on the level of single events must 
necessarily appear as "ad hoc". In contrast, from the viewpoint of a contextual description, 
the introduction of such a model is a necessity. ^^^ 

4. Simulation model 



In our simulation model, every essential component of the laboratory experiment such as 
the source, the Fresnel biprism, and detector array has a counterpart in the algorithm. The 
data is analyzed by counting detection events, just as in the laboratory experiment. ^^^ The 
simulation model is solely based on experimental facts. 

The simulation can be viewed as a message-processing and message-passing process routing 
messengers through a network of units that processes messages. The processing units play the 
role of the components of the laboratory experiment and the network represents the complete 
experimental set-up. We now specify the operation of the basic components of the simulation 
model in full detail. Other components that are specific to a particular interference experiment 
are described together with the presentation of the simulation results. 
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Fig. 4. Diagram of the event-based detector model defined by Eqs. (7) and (8). The dashed line 
indicates the data flow within the processing unit. 



4-1 Messenger 

In our simulation approach, we view each photon as a messenger carrying a message. 
Each messenger has its own internal clock, the hand of which rotates with frequency /. As the 
messenger travels from one position in space to another, the clock encodes the time-of-flight 
t modulo the period 1//. The message, the position of the clock's hand, is most conveniently 
represented by a two-dimensional unit vector e^ = (eo^fc,ei^^.) = (cos 0^, sin 0^), where the 
subscript k > labels the successive messages, (pk = '^'^ftk, and t^ is the time-of-flight of the 
fe-th messenger. The messenger travels with a speed c/n where c denotes the speed of light in 
vacuum and n is the refractive index of the medium in which the messenger moves. 

4-2 Source 

In a simulation model in which the photons are viewed as messengers, the single-photon 
source is trivially realized by creating a messenger and waiting until its message has been 
processed by the detector before creating the next messenger. This ensures that there can be 
no direct information exchange between the messengers, even if one particle has arrived at the 
detector while another particle is at the source or at a detector, implying that our simulation 
model (trivially) satisfies Einstein's criterion of local causality. 

For the double-slit, two-beam interference, and circular slits simulations, messengers leave 
the source at positions generated randomly according to the current distributions Eqs. (1), 
(3), and (5), respectively. The distribution of the angle /3 is chosen to be uniform. When 
messenger k is created, its internal clock time tk is set to zero. 

4.3 Detector 

A single photon detector, such as a photographic plate, consists of many identical detection 
units each having a predefined spatial window in which they can detect photons. Because these 
small detection units are photon detectors themselves we also name them detectors in what 
follows. Here we construct a processing unit that acts as a detector for individual messages. 
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A schematic diagram of the unit is shown in Fig. 4. The first stage consists of a deterministic 
learning machine (DLM) that receives on its input channel the kth message represented by 
the two-dimensional vector e^ = {cos (j))., sin (p^). In its simplest form the DLM contains a 
single two-dimensional internal vector with Euclidean norm less or equal than one. We write 
Pfc = {PQ,k-,Pi,k) to denote the value of this vector after the fcth message has been received. 
Upon receipt of the fcth message the internal vector is updated according to the rule 

Pfc = 7Pfc-i + (1 -7)efc, (7) 

where < 7 < 1 and /c > 0. Update rule Eq. (7) clearly indicates that the first stage 
learns from the incoming messages in a deterministic way and therefore it is given the name 
deterministic learning machine. Obviously, if 7 7^ 0, a machine that operates according to the 
update rule Eq. (7) has memory. 

The second stage of the detector (see Fig. 4) uses the information stored in the internal 
vector to decide whether or not to generate a click (threshold behavior). As a highly simplified 
model for the bistable character of the real photodetector or photographic plate, we let the 
machine generate a binary output signal Sk using the intrinsic threshold function 

5fc = e(pi-rfc), (8) 

where 0(.) is the unit step function and < r^ < 1 is a uniform pseudo-random number. 
Note that in contrast to experiment, in a simulation, we could register both the Sk = ^ 
and S'fc = 1 events such that the number of input messages equals the sum of the S'jt = 
and S'fc = 1 detection events. Since in experiment it cannot be known whether a photon has 
gone undetected, we discard the information about the S'fc = detection events in our future 
analysis. 

The total detector count is defined as 

fc 
N = Y,S,, (9) 

i=i 
where k is the number of messages received. Thus, N counts the number of one's generated by 
the machine. As noted before a detector screen is just a collection of identical detectors and 
is modeled as such. Each detector has a predefined spatial window within which it accepts 
messages. 

In Appendix A we prove that as 7 — t- 1~, the internal vector p^ converges to the average 
of the messages ei, . . . , e^. In general, the parameter 7 controls the precision with which the 
machine defined by Eq. (7) learns the average of the sequence of messages ei,e2,... and 
also controls the pace at which new messages affect the internal state of the DLM (memory 
effect). ^^•' In Appendix B we show how to modify the update rule Eq. (7) such that the 
transient regime of the detector becomes shorter. The transient behavior of the simplest 
and the slightly more complicated detector models may be accessible to real experiments, as 
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explained in Section 6. In appendix B, we also give an alternative for Eq. (8) that does not 
make use of pseudo-random numbers. 

Before we proceed we make a few notes on the memory and threshold behavior of our 
detector simulation models. Although the word memory may give the impression that the 
detector keeps track of all the photons that pass, all the event-based detector models intro- 
duced in this paper have barely enough memory to store the equivalent of one message. Thus, 
these models derive their power, not from storing a lot of data, but from the way they process 
successive messages. Most importantly, the DLMs do not need to keep track of the number 
k of messages that they receive, a number that we cannot assume to be known because in 
real experiments we can only count the clicks of the detector, not the photons that were not 
detected. As shown in Appendix C, the role of the local memory of the detector is similar to 
that of the dielectric function in Maxwell's theory. Our detector models do not incorporate 
a memory fade-out as a function of time. Although this could be an essential feature in time 
frames in which the detectors do not receive photons, we do not consider it to be of importance 
for our present study. 

We also want to emphasize that the presence of a threshold does not cause our detector 
model to operate with less than 100% efficiency. In general, the detection efficiency is defined 
as the overall probability of registering a count if a photon arrives at the detector^'^) Using this 
definition, our event-based detector model simulates an ideal single-photon detector that has 
100% detection efficiency. This can easily be demonstrated by performing the simulation of an 
experiment (which is very different from a double-slit experiment) that measures the detection 
efficiency. ^^^ In such an experiment a point source emitting single photons (messengers) is 
placed far away from a single detector. As all photons that reach the detector have the same 
time-of- flight (to very good approximation), all the messengers that arrive at this detector 
will carry the same message. As a result, the internal vector rapidly converges to one, so that 
the detector clicks every time a photon arrives. Thus, the detection efficiency, as defined for 
real detectors, of our detector model is very close to 100%. Although the detection efficiency 
of the detector model itself is very close to 100%, the ratio of detected to emitted photons 
is much less than one. Note however that, in general, as is well known, a photon detector 
+ electronics is an open system (powered by external electrical sources etc.), hence photon- 
energy conservation within the detector-photon system is not an issue. 

4-4 Discussion 

In our approach, interference appears as a result of processing individual events, but 
definitely not because we have introduced "wave-like" ingredients in a sneaky manner. In our 
corpuscular model, each particle carries its own clock, that is, it carries its own local oscillator. 
This oscillator only serves to mimic the frequency of the individual particle (photon). A the 
particle hits the detector, the detector "observes" the state of the oscillator that is attached to 
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this particular particle and determines its time-of- flight. Note that the idea of introducing the 
time-of-flight does not mean that we obtain interference by summing wave functions a^e"*'^** 
where tk denotes the time-of-flight of the A;th particle. 

There is no communication/interaction between the detectors that make up the detection 
screen, hence there is no wave equation (i.e. no partial differential equation) that enforces 
a relation between the internal variables of these detectors. Likewise, the oscillator that is 
carried by a particle never interacts with an oscillator of another particle, hence the motion of 
these two oscillators is also not governed by a wave equation. Naively, one might imagine the 
oscillators tracing out a wavy pattern in space as they travel from the source to the detector 
screen. However, in our model there is no relation between the times at which the particles 
leave the source, hence it is impossible to characterize all these traces by a field that depends 
on one set of space-time coordinates, as required for a wave theory. 

5. Simulation results 

First, we demonstrate that our event-by-event simulation model reproduces the wave 
mechanical results Eq. (2) of the double-slit experiment. Second, we simulate a two-beam 
interference experiment and show that the simulation data agree with Eq. (4). Third, we 
validate the simulation approach by reproducing the interference patterns for two circular 
sources, see Eq. (6). Finally, we present the results for the simulation of the single-photon 
interference experiment with a Fresnel biprism,^^-' see Fig. 1. The results presented in this 
section have all been obtained using the detector model described in Section 4. Simulation 
data produced by the detector models described in Appendix B are given in Section 6. 

5.1 Double-slit experiment 

As a first example, we consider the two-slit experiment with sources that are slits of width 
a = A (A = 670 nni in all our simulations), separated by a center-to-center distance d = 5A, see 
Fig. 2. In Fig. 5(a), we present the simulation results for a source-detector distance X = 0.05 
mm. When a messenger (photon) travels from the source at (0, y) to the circular detector 
screen with radius X, it updates its own time-of-flight, or equivalently its angle (p. This time- 
of-flight is calculated according to geometrical optics. ^^' More specifically, a messenger leaving 
the source at (0, y) under an angle /? (see Fig. 2) will hit the detector screen at a position 
determined by the angle 9 given by 



sin0 = zcos /? -|- sin/^Y 1 — z^cos^ /3, (10) 

where z = y/X and \z\ < 1. The distance traveled is then given by 



s = Xy/l-2zsmf3 + z^, (11) 

and hence the message is determined by the angle (p = lufsjc where c is the speed of light. 
As the messenger hits a detector, the detector updates its internal vector and decides whether 



11/30 



J. Phys. Soc. Jpn. 



Full Paper 




Fig. 5. Detector counts as a function of the angular (spatial) detector position (y) as obtained from 
event-by-event simulations of the interference experiment shown in Fig. 2 (Fig. 3). The circles 
denote the event-based simulation results produced by the detector model defined in Section 4. 
The dashed lines are the results of wave theory (see Eqs. (2), (4) and (6)). (a) The sources are slits 
of width a = A (A = 670 nm in all our simulations) , separated by a distance d = 5A and the source- 
detector distance X = 0.05 mm, see Fig. 2. The sources emit particles according to the current 
distribution Eq. (1). An interactive program for the double-slit simulation can be downloaded 
from the Wolfram Demonstration Project web site;^^^ (b) The sources 5i and 52, separated by a 
distance d — 8A, emit particles according to a Gaussian current distribution Eq. (3) with variance 
a — X and mean d/2 and —d/2, respectively (see Fig. 3). The source-detector distance X = 0.1 
mm; (c) The two circular sources 5*1 and ^2 of radius a ~ X with centers separated by a distance 
d = 5A emit particles according to the current distribution Eq. (5). The distance between the 
center of the two-source system and the spherical detection screen is X = 0.1 mm. 



to output a zero or a one. 

This process is repeated many times. The initial y-coordinate of the messenger is chosen 
randomly from a uniform distribution on the interval [—d/2 — a/2, —d/2 + a/2] U [+d/2 — 
a/2, +d/2 + a/2]. The angle /3 is a uniform pseudo-random number between —it/2 and 7r/2. 

The markers in Fig. 5(a) show the event-by-event simulation results produced by the 
detector model described in Section 4 with 7 = 0.999. We used a set of thousand detectors 
positioned equidistantly in the interval [—57°, 57"], each of them receiving on average by six 
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thousand photons. The number of chcks generated by the detectors, that is the number of 
so-called detected photons, is approximately 16.10^. Hence, the ratio of detected to emitted 
photons is of the order 0.25, a fairly large number compared to those achieved in laboratory 
experiments with single-photons (see Section 5.4). The result of wave theory, as given by the 
closed- form expression Eq. (2), is represented by the dashed line. Without using any knowledge 
about the solution of a wave equation, the event-based simulation (markers) reproduces the 
results of wave theory. 

According to our mathematical analysis of the performance of the machines (see Appendix 
A), accurate results (relative to the predictions of quantum theory) are to be expected for 
7 close to one only. Taking for instance 7 = 0.99 does not change the qualitative features 
although it changes the number of counts by small amounts (data not shown). 

An interactive Mathematica program of the event-based double-slit simulation which al- 
lows the user to change the model parameters and to verify that the simulation reproduces 
the results of wave theory may be downloaded from the Wolfram Demonstration Project web 
site.29) 

5.2 Two-beam interference experiment 

As a second example we consider the two-beam interference experiment depicted in Fig. 3. 
We assume that the messengers leave either source Si or S2 from a position y that is distributed 
according to a Gaussian distribution with variance a and mean +d/2 or —d/2, respectively. 
Also in this case, the time-of-flight is calculated according to geometrical optics. ^''•' A messen- 
ger leaving the source at (0, y) under an angle /? (see Fig. 3) will hit the detector screen at a 
position {X, y') 

y' = Xtan^ + y, (12) 

the distance traveled is given by s = Xsec/3 and the message is determined by the angle 
(p = 2iTfs/c where c is the speed of light. 

The simulation results for a source-detector distance X = 0.1 mm, for 7 = 0.999 are shown 
in Fig. 5(b). The dashed line is the result of wave theory, see closed form expression Eq. (4). 
Also in this case, the agreement between wave theory and the event-by-event simulation is 
extremely good. 

5.3 Double-slit experiment with circular sources 

As a third example, we consider the double-slit experiment with circular sources, a 
straightforward extension of the two-dimensional double-slit system to three dimensions. As 
shown in Fig. 5(c), there is excellent agreement between the event-by-event simulation and 
the analytical expression Eq. (6). 
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Fig. 6. Schematic diagram of the sim.ulation setup of a single-photon experim.ent with a Fresnel 
biprism. The apex of the Fresnel biprism with summit angle a is positioned at (X',0). In the 
simulation, a line source emits particles from positions drawn from the current distribution Eq. (13) 
with random angles /? chosen uniformly from the interval [—a/2, a/2]. The detectors D positioned 
at {X, y) count the photons. 



5.4 Experiment with a Fresnel biprism 

Finally, we consider the single-photon experiment with a Fresnel biprism. ^^' Figure 6 shows 
the schematic representation of the single-photon interference experiment that we simulate. 
For simplicity, we assume that the source S is located in the Fresnel biprism. Then, the results 
do not depend on the dimensions of the Fresnel biprism. Simulations with a Fresnel biprism 
of finite size yield results that differ quantitatively only (results not shown) . 

Messengers are created at positions drawn randomly from the distribution 

J{x,y) = 6{x)e-y'/^'^\ (13) 

As in all other cases, the time-of-flight of the messenger is calculated according to the rules of 
geometric optics. ^''•' A messenger starting at (0, y) with angle /3 (see Fig. 6) leaves the Fresnel 
biprism at 

X' =F y tan a/2 



x± 



y± 



1 lb tan /3 tan a/2 ' 
y' + X' tan 13 



(14) 



1 lb tan /3 tan a/2 ' 

where the sign has to be chosen such that x± < X' and ^y± > 0, that is such that the path 
of the messenger crosses the Fresnel biprism boundary. Using the fact that the tangential 
component of the velocity is continuous across the Fresnel biprism boundary, ^"^^ we have 



P'± 



iba 



-b arcsin 



nsin(/3 =F 



a. 



(15) 



and we find that the messenger hits the screen at D = (X, {X — x±) tan/3^ -b y±) and that 
the total time traveled is given by 

x± 



n — seep H 

c c 



■ sec j3'j^ . 



(16) 
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Fig. 7. Detector counts as a function of the detector position y of the detector array positioned at 
X for a single photon interference experiment with a Fresnel biprism (see Fig. 6). The Fresnel 
biprism has an index of refraction n — 1.5631 and a summit angle a = 1°. Its apex is positioned at 
(X', 0) with X' — 45 mm. The source emits particles according to a Gaussian current distribution 
with variance a — 0.531 mm and wavelength A = 670 nm.^^^ The circles denote the event-based 
simulation results. The dashed lines denote the numerical results as obtained from wave theory. 
(&) X — X' = 1 mm; (b) X — X' ~\^ mm; (c) X — X' — h^ mm. Thousand detectors where used 
to record the individual events. 



In the simulation, the angle of incidence /3 of the photons is selected randomly from the 
interval [— Q/2,a/2], where a denotes the summit angle of the Fresnel biprism. A collection 
of representative simulation results for 7 = 0.999 is presented in Fig. 7. The dashed lines are 
the numerical results obtained from wave theory by Monte Carlo sampling. Again, we find 
that there is excellent quantitative agreement between the event-by-event simulation data and 
wave theory. Furthermore, the simulation data presented in Fig. 7 is qualitatively very similar 
to the results reported in Ref.^^^ (compare with Fig. 4(d) and Fig. 5(a) (b) of Ref.^^-*). Figure 
4(c) and (4d) of Ref.^^^ are made of approximately 20000 photocounts on the CCD camera, 
while the number of photodetections on the avalanche photodiodes in absence of the CCD 
camera would be 40.10^ during the exposure time of 2000 s. Hence the ratio of detected to 
emitted photons is of the order of 0.0005. This ratio is much smaller than what we observe in 
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our idealized simulation experiment. Namely, each of the thousand detectors making up the 
detection area is hit on average by sixty thousand photons and the number of clicks generated 
by the detectors is approximately 16.10^. Hence, the ratio of detected to emitted photons is 
of the order of 0.026, much larger than the 0.0005 observed in experiment.^^-* 
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Fig. 8. (color online) Detector counts as a function of the angular detector position 6 for the inter- 
ference experiment shown in Fig. 2 which employs only one single-photon detector that is swept 
over the half circle with a fixed angular velocity. The results are obtained from event- by- event 
simulations with six different detector models. The line sources have a width a — X are separated 
by a center-to-center distance d — 3A, and X — 0.05 mm (see Fig. 2). The labels in the figures 
indicate the detector model (algorithms) used. Roman numbers refer to the DLM update rule. 
The letters a and b refer to the pseudo-random and deterministic generation of clicks, respectively. 
la: Eqs. (7) and (8); lb: Eqs. (7) and (24); Ila: Eqs. (22) and (8); lib: Eqs. (22) and (24); Ilia: 
Eqs. (23) and (8); Illb: Eqs. (23) and (24). Stars: Nsweeps ~ f; Crosses: Nsweeps ~ 25; Triangles: 
^sweeps = 50; Diamonds: A^siueeps = 100; Solid lines: Wave theory, see Eq. (2). Other lines are 
guide to the eye only. 
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6. Experimental tests: A proposal 

The simulation models that we propose in this paper make specific predictions that may be 
tested by carefully designed, time-resolved single-photon interference experiments. However, 
not all experiments one can think off are as easy to realize. One of the simplest proposals to test 
the simulation models would be to consider a large number (M) of identical and independent 
two-beam (or double-slit) interference experiments in which only one photon is detected at 
each of the M detection screens. According to quantum theory, summing up the single spots of 
the M detection screens gives the same interference pattern as if one would conduct one two- 
beam interference experiment with M photons being detected on the same detection screen 
(all under the assumption that every time a photon is emitted and that all emitted photons 
are detected). For this experiment, the simulation models that we have introduced do not 
yield an interference pattern, as is clear from their description. Thus, at least in principle, 
this experiment should be able to refute the corpuscular model. Note that in the absence 
of any experimental evidence and bearing in mind that quantum theory has nothing to say 
about individual events, ^^^ it is only a hypothesis that the experiment with finite M will yield 
results that agree with quantum theory. Whether this hypothesis is actually true remains to be 
demonstrated by an experiment. Unfortunately, in practice, this experiment may be difficult 
to realize, the central question being how large M should be before one observes a pattern 
that resembles the one predicted by wave theory. A rough estimate, based on experiments 
with electrons^^' suggests that M > 50000, a number which makes this proposal very hard 
to realize in practice. Therefore, we propose another experiment that may be realizable with 
present-day technology. 

As explained earlier, if our simulation models operate in the stationary-state regime, they 
reproduce the wave theoretical results. Therefore, to falsify our event-based models the single- 
photon experiment should be designed such that it is sensitive to the transient behavior 
of the whole setup. In other words, the experiment should operate on a time scale that is 
sufficiently short to prevent the DLM in our detector models to reach the stationary state. 
For a fair comparison between experiment and our simulation models, it is essential that the 
experimenter does not discard data that is recorded during the "calibration" or "warm-up" 
stage because this data may contain valuable information about the transient behavior of the 
experimental setup. 

In this section, we use our simulation approach to make predictions of laboratory exper- 
iments that may be realizable. Consider again the double-slit experiment depicted in Fig. 2 
but instead of having many detectors at different angles 9, we use only one detector placed 
on a goniometer. The idea is to keep the total exposure time constant while the detector is 
swept back-and- forth over (part of) the half-circle (see Fig. 2). In our simulation models, the 
recorded interference pattern will then depend on the angular velocity of the detector. For 
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velocities that are sufficiently small to allow the DLM to reach the stationary state, the inter- 
ference pattern obtained agrees with the one predicted by wave theory. On the other hand, if 
the detector position changes rapidly, the DLM may not receive enough events to accurately 
reproduce the wave mechanical result. Therefore, if we keep the total exposure time constant 
and perform a set of experiments for several choices of the sweep velocity, our simulation 
models predict that the interference patterns will change and that these changes reflect the 
internal dynamics of the detector model used. 

The procedure that we propose is the following. First, we fix the angle 56 by which the 
detector position will be moved. For simplicity, we assume that the aperture of the detector 
is equal to 59. Then, we fix the total number of events N total which, on average, will arrive 
within each arc of angle 50. Finally, we select the number of times Ngweeps that the detector 
will be swept back-and-forth over the half circle. 

In the simulation, the internal variables of the detector models are initialized once. The 
simulation results presented in this section have been obtained using 59 = 1°, N total = 10^, 
^sweeps = 1,25,50,100, 7 = 0.999, and for the modified detector models introduced in Ap- 
pendix B, K = 0.9, wq = 0.9 and u = 0.99. In all figures, the theoretical result Eq. (2) is 
rescaled to fit to the maximum of the simulation data at the smallest sweep velocity and, in 
the case of Ilia and Illb, also shifted to account for the non-zero bias. 

As explained in Appendix B, the simple detector model introduced in Section 4 with 
the DLM defined by Eq. (7) may require a significant amount (order of thousands) of input 
events to reach the stationary state. Hence, if we move the detector before the DLM reaches 
its stationary state, this detector model may not produce results that agree with wave theory. 
This expectation is confirmed by the results shown in Figs. 8(Ia) and (lb). If the detector 
moves slowly {N sweeps = 1)) the event-based simulation data are in concert with wave theory, 
as is clear from the comparison of the stars and the solid lines in Figs. 8(Ia) and (lb). From 
Figs. 8(Ia) and (lb) it is also clear that increasing the number of sweeps to Ngweep = 25 (recall 
that the total amount of events corresponding to the total exposure time in the experiment is 
fixed) leads to a reduction of the visibility of the fringes. If we increase the number of sweeps 
to Nsweep = 50, the detector model fails qualitatively. 

Thus, an experiment that uses a moving detector might be able to rule out event-based 
models la and lb as candidate descriptions of the single-photon interferences. However, this 
does not yet imply that our approach as such should be abandoned: It may be that the 
detector model is too simple. Therefore, it is of interest to explore to what extent the results 
depend on the particular algorithms used. 

It is not difficult to modify the DLM defined by Eq. (7) such that the convergence to the 
stationary state is much faster or that the response to changes in the input data is faster. In 
Appendix B, we give the details of two of such variants. 
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DLM II is constructed such that its stationary state behavior is the same as that of 
the simple DLM (Eq. (7)), hence the detector model using this DLM reproduces the results 
of wave theory if we employ an array of detectors or move the single detector very slowly. 
From Figs. 8(IIa) and (lib), we may conclude that this model is an improvement over the 
simple model in that it still shows interference fringes at a sweeping rate of Ngweep = 50. For 
^sweeps = 100, the detector receives approximately N total / iX'^^N sweeps / ^G) ^ 55 events before 
it moves to the next position. With this small amount of input events, DLM II does not reach 
the stationary state (see also Fig. 9). 

DLM III is a little different than DLM II: It is sensitive to differences between the internal 
state and the input message. As Figs. 8(IIIa) and (Illb) show, these detector models produce 
output signals that are insensitive to the speed at which the detector moves but this comes at 
the price of a nonzero bias which is, within statistical fluctuations, independent of the detector 
position or velocity. Subtracting this bias, all the data fit the theoretical curve very well. 

Summarizing: For experiments that use detectors that have fixed positions, our event- 
based models for the detector yield results that cannot be distinguished from those of wave 
theory. However, our simulation models for single-photon two-beam interference show features 
that may be tested experimentally by measuring the intensity as a function of the speed of a 
moving detector. We have proposed and analyzed a realizable, time-resolved experiment that 
directly probes the dynamics of our detector models and predicted the outcome of such future 
experiments. 

7. Conclusion 

We have demonstrated that it is possible to give a corpuscular description for single- 
photon interference experiments with a double-slit, two beams, and with a Fresnel biprism. 
Our event-by-event simulation model 

• does not require any knowledge about the solution of a wave equation, 

• reproduces the results from wave theory, 

• satisfies Einstein's criterion of local causality, 

• provides a simple, logically consistent, particle-based description of interference. 

We do not exclude that there are other event-by-event algorithms that reproduce the 
interference patterns of wave theory. For instance, in the case of the single-electron experiment 
with the biprism, '^^' it may suffice to have an adaptive machine handle the electron-biprism 
interaction without having adaptive machines modeling the detectors. We leave this topic for 
future research. 

We hope that our simulation results will stimulate the design of new time-resolved single- 
photon experiments to test our corpuscular model for interference. In Section 6, we proposed 
such an experiment and also predicted the outcome if our simulation model captures the 
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essence of the event-based processes. Note however that the models we have employed are not 
unique, as shown explicitly in Section 4. This leaves some freedom to adapt the simulation 
models to the actual experiments that will be performed. 

Finally, it may be of interest to mention that our approach opens a route for incorporating 
interference phenomena into ray-tracing software. In our simulation method, each messenger 
simply follows one of the rays through the medium, updating the message (corresponding to 
the phase information) as it travels along. Therefore, for applications where the solution of 
the Maxwell equations is prohibitive, the combination of our technique and ray tracing may 
be a viable alternative. 
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Appendix A 

We demonstrate that as 7 — t- 1~ the internal vector p^ in Eq. (7) converges to the average 
of the messages ei, e2, . . .. 

Let ||x|| denote the Euclidean norm of the vector x. Then, as < 7 < 1, ||e,fc|| = 1 for all 
k > 0, and ||po|| = 1 it follows immediately from Eq. (25) that ||pfc|| < 1 for all k > 0, hence 
liuik-^QoPk exists. To determine p = lim^-^ooPfci we have to make assumptions about the 
properties of the sequence {ei, e2, . . .}. For instance, if the sequence {ei, e2, . . .} is generated 
by a stochastic process with mean (ej+i) = e for j = 0, . . . ,k — 1, then it is easy to show that 
p = e. Thus, in this case, the machine defined by the rule Eq. (7) learns the average e by 
updating its internal vector for each message it receives. 

In practice, only finite sequences {ei, e2, . . . , e^^} are available. In this case, we can esti- 
mate the limiting value by assuming that the sequence repeats itself, an assumption that is 
common in Fourier analysis and signal processing in general. ^^•' From Eq. (25), we have 

riK-l 

PnK = T^'Pln-Di^ + (1 - 7) E r^'-'-'ej + l 

j={n-l)K 

K-l 

= 7^P{n-l)K + (1 - 7) X] 7'^"^"^ej-+l+(n-l)K 

i=o 
= 7S(n-i)i^ + (l-7)fe, (17) 



where 



K-l 

fi^ = E^''''"Vi> (18) 

i=o 
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and n > 0. From Eq. (17) we find 



-1 



^nK = t'^^PO + (1 - l) \ ^\ iK, (19) 



and hence 

lim pnK = -, \ Y. T^'-^-'e.+i, (20) 

n— )-oo 1 — O" — 

' J=0 

such that 

1 ^'^ 
Hm hm pnK = 17 V] ^j+i- (21) 

From Eq. (21), we conclude that as 7 ^ 1~ the internal vector p^ converges to the average 
of the messages ei, . . . , ex- In general, the parameter 7 controls the precision with which the 
machine defined by Eq. (7) learns the average of a sequence of messages and also controls the 
pace at which new messages affect the internal state of the learning machine.^ ' 

Appendix B 

Without performing any simulation, we can already see from Eq. (7) that the simple 
machine may not perform very well in some cases. Suppose that po = and that e^ = e for 
all k. Then, from Eq. (7) it follows that p^ = (1 — 7'^)e such that ||pfc — e|| = 7^. Although 
the latter equation shows that the convergence of p^ to the input vector e is exponentially 
fast, for 7 very close to one, in practice, it may take quite a number of events to reach the 
stationary state. 

In this Appendix, we describe two modifications of the algorithm Eq. (7) of the first stage 
(DLM) and one alternative for the algorithm Eq. (8) of the second stage. The modifications 
of the first stage reduce the amount of events required for the detector model to reach the 
stationary regime. The alternative for the second stage eliminates the need for a pseudo- 
random number generator. 

It is not difficult to modify the machine such that its asymptotic behavior remains the 
same while improving, significantly, the speed with which it learns from the input ej^.. A simple, 
but by no means unique, modification is to add one memory element to store one variable, 
denoted by Wk, which keeps track of the differences between p^ and pfc-i. For k > 0, these 
variables are updated according to the rule 

Pfc = Mfc-iPfe-i + (1 - ^fc-i)efc, 

Wk = KWk-l + {l-K) , (22) 

where < k < 1 is another control parameter and < wo < I- Although the variable fik 
is redundant, we wrote Eq. (22) such that it is obvious that it is an extension of Eq. (7). In 
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Fig. 9. (color online) The square of the length of the internal vector p^ as a function of the number 
of received events k for three different input messages e^ = (r^, , (1 — r^)^/^) (top lines), e^ = 
(cos Trrfc , sin Trrfc ) (middle lines), and e^ = (cos27rrfc, sin27rrfc) (bottom lines) where the < r^ < 1 
are uniform pseudo-random numbers. Dashed lines: Model Eq. (7). Solid lines: Model Eq. (22) 
with K = 0.9. The inset shows the short-time response of models Eq. (7) and Eq. (22) in more 
detail. In all cases po — (1, 0), wq — 0.9 and 7 ~ 0.999. 



essence, instead of keeping 7 fixed in the rule to update p^ (see Eq. (7)), in Eq. (22), the value 
of /ifc in the rule to update p^ is made variable. This flexibility is then exploited through the 
first and last rule in Eq. (22). The last rule defines a machine that learns the distance between 
Pfc and pfc-i, the learning speed being controlled by n. The basic idea is that if this distance 
is large (say close to but less than 2), the last rule will drive Wk to one such that fi^ is small 
and the change of p^ may be large. In the opposite situation, the last rule will force Wk to zero 
and pfc will change by small amounts (assuming 7 is close to but less than one). As yU^ < 7, 
the asymptotic behavior of the machine defined by the rule Eq. (22) is easily shown to be the 
same as that of the simple version in which we keep /U^ = 7. Thus, although equations that 
govern the dynamics of the machine Eq. (22) are nonlinear (in the p's), asymptotically the 
dynamics is governed by the linear equation Eq. (7). 

It is not easy to study the transient behavior of the classical, dynamical systems defined 
by Eqs. (7) and (22) by analytical methods but it is almost trivial to simulate these models 
on a computer. In Fig. 9, we show some representative simulation results to illustrate that 
the slightly more complicated machine Eq. (22) performs significantly better than the simple 
machine Eq. (7) with respect to the number of events it takes for the machine to reach 
the stationary state. Roughly speaking, after about 60 events, machine Eq. (22) has learned 
enough to reproduce the correct averages. As expected on theoretical grounds, both machines 
converge to the same stationary state. 
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A minor modification of algorithm Eq. (22) yields the DLM defined by 

Pfc = fJ'k-iPk-i + (1 - Aifc-i)efc, 

/I \ Pfc ^k /oo\ 

Wk = KWk-l + {l-K) . (23) 

Note that the only change is in the third rule where we replaced Pfc_i by e^. This replacement 
causes the machine to respond very fast to changes in the sequence of input messages {e^} 
but, at the same time, also leads to a reduction of the average value of Hk-i which in turn, 
will cause the detector model to produce a nonzero signal, independent of the input messages 
(see Figs. 8(IIIa) and (Illb)). 

As an alternative to the pseudo-random "click generator" Eq. (8), we may generate the 
clicks by means of a very simple DLM^^' containing a single internal variable < z^ < 1 that 
is updated according to 

^k = \ , 

[ 1 otherwise 

Zk = vzk-i + (1 - v)Sk. (24) 

Here, the parameter {) < u <1 plays the same role as 7 in Eq. (7). The non negative number 
p| is the input message for the DLM. The dynamics of the system defined by Eq. (24) is very 
different from that of Eq. (7).^^^ Elsewhere, we have shown that for a fixed input message p^, 
the machine defined by Eq. (24) generates a binary sequence (the Sks) such that in the long 
run the ratio of the number of ones relative to the total number of events is equal to the time 
average of p^-^^^ Thus, the machine defined by Eq. (24) produces clicks with a rate that is 
determined by p^. 

Appendix C: Relation between simulation model and wave mechanics 

The simulation results presented in Section 5 demonstrate that the event-based model is 
capable of reproducing the results of wave theory without making recourse to the solution of 
the wave equation or even a single concept of wave theory. As there seems to be a general 
consensus that such models are not supposed to exist, it is of interest to show that for the 
problems that we deal with in this paper, the event-based model contains the description that 
derives from Maxwell's equations. 

Our demonstration consists of two steps. First we relate the variables of the event-based 
model to those of classical electrodynamics. Second, in analogy with the derivation of the 
diffusion equation from the discrete random walk model, we show how our event-based model 
leads to the Debye model for the interaction between material and electric field. Other models 
such as the Drude or Lorentz model can be derived in a similar manner but to keep the 
presentation concise, these derivations are relegated to a future paper. 
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Table I. Correspondence between Maxwell's theory and the particle-based, event-by-event simulation 
model. For simplicity of presentation, we consider the case of a linearly polarized wave only. 

Classical electrodynamics Event-based simulation model 



Description 



wave 



particle 



Properties 



oscillator frequency / 

direction q 

propagation time t 

phase velocity c 



oscillator frequency / 

direction q 

time-of-flight t^ 

velocity c 



Message 


E = Eq cos(wt — q • r -I- 99) 


ek = (cos27r/ifc,sin27r/tfc) 


Material 


Polarization P(t) 


Internal vector p^ 


Interaction 
with material 


P{t) = f;^x{u)Bit-u)du 


Pfe = 7Pfe-i + (1 - 7)efc 



As is evident from Table I, the messenger can be viewed as the event-based equivalent of 
a classical, linearly polarized electromagnetic wave with frequency /: The message e^ corre- 
sponds to a plane wave with wave vector q (g = 27r//c). The time-of-flight t^ corresponds to 
the phase of the electric field. Adding another clock to the messenger suffices to model the sec- 
ond electric field component orthogonal to the first one, and hence the fully polarized plane 
wave.^'^^ For the systems studied in the present paper including this extra feature, namely 
the equivalent of the polarization of the wave, is not necessary and therefore we confine the 
discussion to messages that are represented by two-dimensional unit vectors. 

The internal vector p^ plays the role of the polarization vector P(t) of the detector ma- 
terial. Indeed, comparing the formal solution of Eq. (7) 



fe-i 

Pfc = T^'Po + (1 - 7) X] ^^^k~j, 
i=o 



with the constitutive equation 



P(t) 



x{u)'E{t — u)du, 



(25) 



(26) 



in Maxwell's theory,^''-' it is clear that both equations have the same mathematical structure: 
The left hand sides are convolutions of the incoming (applied) message (field) with memory 
kernel 7-' (x('w)) (iii applications, we may assume that the initial value po = 0). Thus, the 
DLM is a simple model for the interaction of the individual photons with the material of the 
detector. The time-of-flight, corresponding to the phase of the electric held, is used to update 
the internal vector which corresponds to the polarization vector of the material. 
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Table II. Analogy between the derivation of the diffusion equation from the random walk model and 
the derivation of one of the constitutive equations in Maxwell's theory from the discrete model 
Eq. (7) proposed in this paper. The assumptions that the limiting values D ~ lim^^o linir-^o 5"^ /"^t 
and r = lim^^i- limT-_^o(l ^ 1)1^1 are nonzero and finite are essential to obtain a well-defined 
continuum approximation of the discrete update rules. 



Random walk 



Detector model 



Update rule 

Length scale: 5 
Time scale: r 

Small T 
Small b 
Small b and r 
lim^^o linir^o 
lim^^i- limr^o 

Equation 
Fourier space 



P/,fc+i = \{vi+\,k +Pl-i,k) 
Pl,k =p{lb,kT) =p{x,t) 



Pfc = 7Pfe-i + (1 - 7)efe 
Pfc = p(A;r) = p(t) 
efc = e(A;r) = e(t) 



p,,fc+i = p{x, t) + r%^ + 0(t2) p,_i = Pit) - T^ + 0(r2) 



Pl±l,k =P{x,t) ±b—g^ 

dp(x,t) ^ 5^ d^p{x,t) 



at 

dp(x,t) g2 d^p(x,t) 



+ 



2 dx'-^ 



dt 



2t dx 



+ 0{b^) 

dpjt) 
dt 



1-7 



Pit) + i^e(t) 



(f ^ D, 0<D <oo 



dp{x,t) j-^d'^p{x,t) 

dt dxP 

p(tj) = riiuj + r)-'^ei 



4 

Diffusion equation 



^^F, 0<F<oo 
^ = -rp(t) + Fe(t) 



Picj) = x(^)E(l^) 
Constitutive equation 



Next, we show that this analogy can be carried much further by mimicking the derivation 
that relates the discrete random walk on a line to the one-dimensional diffusion equation. ^^' 
The essential steps for both the random walk and our event-based detector model are sum- 
marized in Table II. Both models describe a process that proceeds in discrete time steps r. 
The random walk model is formulated on a lattice with mesh size b. In the case of the random 
walk, we let the time step r and mesh size b go to zero. In the event-based model we let the 
time step r, that is the time between the arrival of successive messages, approach zero and let 
7 approach one. For both models, we demand that the resulting continuum equations make 
sense. This enforces relations between r and b'^ and between r and 7, as shown in Table II. 
Then, the former relation yields an explicit expression of the diffusion coefficient D = b'^ /2t 
in terms of the length and time scale of the discrete random walk model. Likewise, the latter 
leads to the Debye model for a dispersive medium^^) and gives an explicit expression for the 
relaxation time 1/F = r7/(l — 7) in terms of the parameters of the event-based model. 

As Table II shows, under certain conditions, the discrete models can be approximated by 
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continuum equations that describe the coarse-grained (in space-time for the random walkers 
and in time for the event-based model) behavior but the discrete models provide a description 
with details that can never be extracted from the corresponding continuum equations. Of 
course, the ultimate justification of the event-based model is that, as shown in Section 5, 
it can reproduce the results of wave theory. Appendix D gives a further justification of our 
approach from a computational point of view. 

Appendix D: Computational point of ^riew 

There is a general consensus that unless we first solve the wave equation and use this 
solution as the probability distribution for generating events, there are very fundamental, 
apparently unsurmountable, problems to derive from a wave mechanical description a process 
that produces the events that are observed in experiment. ^^^ The arguments used are rather 
abstract and general^^^ and to understand the subtilities that are involved it may help to 
address this issue from a computational point of view. 

For phenomena that cannot (yet) be described by a deductive theory, it is common prac- 
tice to use probabilistic models. Although Kolmogorov's probability theory provides a rigorous 
framework to formulate such models, there are ample examples that illustrate how easy it is 
to make plausible assumptions that create all kinds of paradoxes, also for every-day prob- 
lems. ^^'^'^^'"^^^ Subtle mistakes such as dropping (some of the essential) conditions, like in 
the discussion of the double-slit experiment, ^^'^^' mixing up the meaning of physical and sta- 
tistical independence or changing one probability space for another during the course of an 
argument, can give rise to all kinds of paradoxes.^^'^^'^^"^^^ For instance, Feynman used the 
double-slit experiment as an example to argue that "far more fundamental was the discovery 
that in nature the laws of combining probabilities were not those of the classical probabil- 
ity theory of Laplace" ,^'^) but this statement has been shown to result from an erroneous 
application of probability theory.^^'^^'^^^ 

By construction, if we use a digital computer to produce numbers as we do in this paper, 
we stay in the domain of elementary arithmetic and we do not have to worry about the 
subtleties of Kolmogorov's probability theory. 

Instead of discussing the apparently unsurmountable problem in its full generality, which 
we could, it is more instructive to examine in detail the simple, concrete example of the 
double-slit model depicted in Fig. 2. According to Maxwell's theory, in the Fraunhofer regime 
the light intensity at the detector on a circular screen is given by^'^) 

(27) 



m 



m 



'— oo 

'sinS^y ^qdsme 
r^ cos^-^ , (28) 

qa sin / 9 
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where p{y') = [0(a — \y' — d/2\) + 0(a — \y' + d/2|)]/2a is the normahzed density distribution 
for the coordinate y' . 

First, starting from the explicit expression Eq. (28) for the density I (9)/ 1(0), it is trivial 
to construct an algorithm that generates events according to this density. Indeed, let us define 

S,{e) = Q{m-r,m), (29) 

where < r j < 1 denotes a uniform pseudo-random number. Then, the number of clicks of 
the detector at angular position 9 is given by 

k 

i=i 
and for sufficiently large k, we have Nii.{9) — > I{9)/I{0) with probability one. This completes 
the construction of the event-based algorithm based on the knowledge of I (9)/ 1(0). Obviously, 
this algorithm is built on the knowledge of the explicit solution 1(9) of the wave problem. 
The events generated by this algorithm build up the interference pattern one-by-one and 
can be identified with the clicks of the detectors. This is as far as the quantum theoretical 
description goes in making contact to the experimental observations: It provides a prescription 
to calculate the probability density to observe a click on a detector. It is quite common to 
postulate that there does not exist a description that goes beyond the specification of the 
probability, excluding that no further advance in a deeper understanding of the process that 
produces the events can be made. 

Disregarding this postulate, we may wonder what happens if we take one step back 
and assume that we only know about expression Eq. (27) in terms of the wave amplitudes 
exp{iqy' sin 9) and density p{y')- Then, the obvious thing to do is to compute the integral in 
Eq. (27) numerically. Without loss of generality, we may write 

^(^) = ]^E^^''''"'' (31) 

^ ' y'&s 
where the summation is over all y' of the set S accounting for the density p{y') and N{S) is 
the normalization factor. By definition of the integral, if the number of elements of the set S 
goes to infinity, we have \A{9)\^ -^ I (9)/ 1(0). 

Although the numerical calculation of the amplitude A{9) is straightforward, there obvi- 
ously is no relation between the points y' of the set S and the number of clicks of the detector 
at 9. In fact, the essence of quantum theory is that there is only a relation between |^(^)| and 
the number of clicks but to know A{9) requires that we first generate (a lot of) pseudo-events 
y'. Obviously, these pseudo-events y' cannot have an interpretation in terms of observed clicks. 

The conclusion therefore is that the description in terms of individual waves (Eq. (27)) 
does not contain the ingredients, not even conceptual, to define a process that generates the 
clicks of the detectors that we observe. Therefore, from a computational perspective, it is futile 
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to try inventing an event-based, particle-like process based on the wave mechanical expression 
for the intensity in terms of sums over amplitudes. 

One may take the position that it is fundamentally impossible to go beyond an event-level 
description based on the knowledge of I (9)/ 1(0) but by postulating this to be true, one simply 
postulates that it is impossible to make any advance in a deeper understanding of event-based 
phenomena. As we have shown by this and many earlier papers, there is no rational argument 
that supports this postulate other than that it is what we have been taught in physics courses. 

Having shown that our event-by-event simulation model reproduces the results of wave 
theory without resorting to a description in terms of waves, we now explain why, from a 
computational point of view, we consider this to be an accomplishment and why our approach 
works. 

The crux of our approach is that we do not start from expression Eq. (27) but construct 
a discrete event process that converges to Eq. (27) while generating events that directly 
correspond to the observed events. During the initial phase, this process may generate events 
that are accidental but once the process has reached its stationary state, the events appear 
with frequencies that corresponds to those predicted by wave theory. 

To understand the idea behind our approach, it may be helpful to draw an analogy with 
the well-known Metropolis Monte Carlo (MMC) method for solving statistical mechanical 
problems. ^^'^^^ The MMC method generates states S, events in our terminology, with a prob- 
ability density ^^' ^^^ 

^-E{S)/kBT 

P^^^l = J2 e--^('S')/fcBT ' (^^) 

where E{S) denotes the energy of the state S, ks is Boltzmann's constant and T is the 
temperature. At first sight, sampling from Eq. (32) is impossible because in all but a few 
nontrivial cases for which the partition function ^^ g-E(S)/kBT jg known, we do not know 
the denominator. The MMC method solves this problem by constructing a Markov chain 
that generates a sequence of events S such that asymptotically these events are distributed 
according to the (unknown) probability density Eq. (32). '^^'^^•' 

The analogy with our approach is the following. Although very different in all details, our 
event-based method uses a deterministic process (implemented as a DLM, see Eq. (7)) of which 
the sampling distribution converges to the unknown probability distribution I (6)/ 1(0). The 
one-to-one correspondence between the objects in the corpuscular, event-based description 
and those in Maxwell's theory (see Section 8) ensures that in the long run, the event-based 
detector model generates clicks with frequencies that correspond to those of the unknown 
probability distribution I{9)/I{0). 
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